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Abstract 

We proposed an evolving network model constituted by the same nodes but different edges. The 
competition between nodes and different links were introduced. Scale free properties have been 
found in this model by continuum theory. Different network topologies can be generated by some 
tunable parameters. Simulation results consolidate the prediction. 
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I. INTRODUCTION 



For the past few years, the interesting properties of complex network, which can be found 
in many fields, such as biology, sociology, physics, and so on, have attracted a number of 
people to explore. Great contribution including the idea and analysis ofrandoar net^i], 
small world networks such as the WS network has been made to the modeling; 



process 



However, real world networks such as Internet, Movie actor and science 
collaboration graph(see 3, [l9( all show a power-law degree distribution property which can 
not be explained by the two models mentioned above. A scale free mo delQQ that satisfies 
a lot of real systems and led many subsequent researchers to further studying has been found 
by Barabasi and Albert in 1999. On the other hand, this model suffers a drawback that the 
exponent of the power law is always fixed, while in real networks it varies. Further study 



has made it explicit that preferential attachment plays a ke y ro 
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In 2000, an exact 



evolving network shows scaling property as expected 
solution to the degree distribution of BA model was obtained by Dorogovtsev, Mendes and 
SamukhinJ^. In their paper, different initial attractiveness is used, to change the exponent 
of the power-law curve. 

Network models are applied to the description of the complex human interactions, known 
as social networks. Enlightened by the fact that there are complicated relationships between 
the elements of social networks, researchers brought multiple-type vertex rather than single 
one, into network n»d*0. Weighted edges as represented by kinds of socia. re.ationsh.ps 
are also introduced in this area[28j. Both models show scaling properties. It is challenging 
to find more fundamental network structures. 

Similarly, our model can be understood from another perspective. One way, but not the 
exclusive one, to understand our concept is to consider the Internet as a network, which 
comprised by online websites as nodes. We set the virtual connections, e.g. the transporta- 
tion forms of information between different websites rather than physical wires and cables, 
as edges between the nodes. It is clear that there are various kinds of communication ways 
between two websites( HTTP, FTP ). So to speak, different kinds of links need to be applied 
to growing networks, where competitions exist between different links as well as the nodes. 
It is attractive to discover the topologies of the whole network and different sub-networks 
containing only one kind of the two edges, and to consider the relation between different 
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edges. In the model we construct, the whole network is somewhat like the partial prefer- 
ence BA model. The significant difference is that we add some tunable parameters into 
the preferential attachment, which makes it possible to alter the exponent of the power-law 
degree distribution and the shift constant of degree distribution function. Our work begins 
with the discussion of a scale free network growing like the BA model with a change of 
the preferential attachment. This is a special example of the growing networks put forward 
by Dorogovstev and Mendes in 2000 26j] in which they add a constant in the preferential 
attachment term. In our main work, the tunable parameters in the model play different 
roles in different edge preferential attachments, and are regarded as a weight to discriminate 
different attractiveness of different links in the evolving network. The research produces 
some interesting results: the whole network and the two sub-networks all evolve just like the 
ordinary scale free network, while the two kinds of links connected to a certain node always 
differ. Besides, the degree distributions of the two sub-networks are unequal. 

II. MODEL AND THEORETICAL APPROACH 

At the beginning, we revise the BA model only a bit and get different expressions of 
the degree distribution function. Before introducing our findings, it is necessary to take a 



look at the work by Dorogovstev and Mendes in 2000 [26[. Their contribution mainly lies 
with the modulation of the original BA model and the way they solved the problem, well 
known as the 'Master Equation Approach'. In their paper, the preferential attachment in 
BA model is no longer fixed; instead, changeable initial attractiveness, which is given to each 
site, together with incoming links of the site, determines the probability whether a new link 
will point to this site. Under long-time limitation, the exponent of the degree distribution 
varies from 2 to oo , depending on the initial attractiveness. In this paper we do not intend 
to discuss too many complicated forms of preferential attachments deriving from the BA 
model. We barely make a simple and clear revision, and attain concise and approximately 
precise results, with the application of continuum theory. 

The preferential attachment in BA model is changed in this way: the possibility of the 
degree increase of a site is proportional to n(^«) = m ^fc + +/) * s ^ e degree of site i, while 
/ is a given constant). This model may be considered as a special example of the model 
generalized by Dorogovestev and Mendes if one takes in their model as m + f in our 
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case. By making use of master equation to solve this problem, the degree distribution was 
found to have the exponent — (2 + a) where a = A linear shift ma = is found in 
the distribution function 

P{q) « 1 + a) 11 \ I(g + ma - (2+a) 1 

1 [ma) 

We here apply continuum theory, rather than master equation, to solve the problem we 
raise. It is not difficult to draw the following equation: 

where m is the number of links every time the network increases. The solution is quite 
simple: as t approaches infinity, the degree distribution follows this equation 

P(k) = (2 + L)( m + f) 2+ ™(f + k)-^ (3) 
m 

Compare the degree functions ((TJ and (|3*]l. taking the power term into account, and we find 
that the two results are the same in essence. In our derivation the tunable parameter is /, 
while in D&M's it's the initial attractiveness . By transforming A^ to m+f or inversely, 
the two functions show the same property. This transformation is reasonable. If we adopt 
D&M's way[26] in our model, A^ should be the initial degree a site has when it enters the 
network, which is exactly m + f. In D&M's paper [26], the exponent 7 = 2 + a varies from 2 
to 00 while in our model, 7 = 3 + ^. / can be taken from — m to 00. Therefore, the range 
is the same. The advantage by using continuum theory is that no special function appears, 
which makes it relatively easy to solve the two-link networks we conceive in the following 
content. Another difference between our model and D&M's model is that we consider the 
network as an undirected network while they consider the network as directed. 

The following is our new proposal based on the idea of a new kind of network. This 
network, which contains two kinds of edges, is very different from the networks we have 
studied till now. If we consider the different links as identical, the model should come back 
to partial preferential attachment BA case. To reach this point briefly, we give different 
linear preferential attachments to two kinds of edges respectively. Of course, modulatory 
parameters are added to make our model changeable to fit the real networks. 

The model: At each time interval, a new node adds to the network and connects to the 
old nodes with m edges. We divide the m edges into two types: one is called X edge (or X 



link) and the other Y edge (or Y link). The probability that node i in the original network 
makes a connection to the new added node follows such preferential attachments: 

dx t = Xj + yj + / + m , , 

dt m 2£x 3 + Vj + f) U 

dm = Xi + Vi + f - 9Vi /,x 
dt m 2(Ex, + Vj + f) U 

where Xi represents the number of X-edges node i has connected at time t; i/i represents the 
number of F-edges. 

Here, we give linear asymmetric expressions in X and Y edge aiming to represent the 
linear relations in two networks. However, we do not assume this to be exactly the reflection 
of real world networks, it is a step forward in the direction of finding out complex relations 
between networks that share same nodes but different links. The parameter g here plays 
the key role in determining the relation between X and Y network: When g is very small, 
it can be neglected so that the two networks evolve in the same way- they just behave as 
the same network. When g approaches infinity, say, some relatively large number, the above 
expressions show that new edges will always be added as X edges while the number of Y 
edges could only remain constant instead of increasing. In general, g is used to determine 
the difference of the increasing rates of the two different edges. 

The number of total edges node i connects satisfies the following equation: 

d(xi + yi) Xi + yi + f 

We have studied this kind of network before, Xj + Hi = ki and degree distribution follows 
equation (jSJ). 

Firstly, we discuss Y sub-network. Because we add m edges at every time step, so 

Y t {x i + y i + f) = (2m + f)t (7) 
According to continuum theory, equation (J3J became 

m 

dyi (m + f)(j-) w ~ 9Vi 



m 



Its solution is 



dt 2(2m + f)t 



771 + f , t . m I rng ■. 

/ Ism i n J- J- U( 2m+ n/ 



(8) 



l/i = J -(-)2m + f _|_ C0nSt ■ t l 2(2m + /)^ (9) 

2 + g ti 



( mg \ 

When t approaches infinity, t ^ 2 ( 2m +/) ; approaches and Xji becomes 

m + f.t rn 
y 2 + <?V 

Then we obtain 



(10) 



P y { Vi <y) = P{ ti > y-^C^r^t)) 

z + g 

1 --(2+^ m+ A (2+;D 



For large t ,P(yi) reads 

P a(s) .^|^ K ( 2 + l) ( ^ + V«^ (12) 

Up to now, we have got the mathematical expression of degree distribution of Y sub-network, 
i.e. equation (fT2|) . The degree distribution is power law with the same exponent as that 
of total degree distribution but without shift /. Similarly, we proceed in X sub- network. 
According to equation (J3J) and fTOj) . we get 

%i ">i Vi 

1 —I— (1 T, r m \ l rnq \ 

= —^-(m + /)(-) ( w) - / -Const -r ( 2TOT)) (13) 

2 + 5 ti 



Under long-time limitation, we neglect Const ■ t 2 ( 2m +fl , and obtain 



l±l {m + f)i h(^-f (M) 
z -+- g tj 

The same as what we dealt with Y sub-network, the degree distribution of X sub-network 
becomes 

P*{x) = (2 + Q [i±|(m + f)f + ^\x + /)-( 3+ ^) (15) 

It is a power-law distribution with shift /. The preferential attachments (jlj) and (|SJ) are not 
the same. Therefore, the number of x edges does not equal to the number of y edges. From 
the calculations above, we get the ratios of total degree to x degree and y degree 

P x+y (k-f):P x (k-f):P y (k) 

= ! +i±* r™ a j_v + * (16) 



2 + gJ \2 + g 

the subscript x + y in P x+y (k) means the degree in total network 



III. NUMERICAL SIMULATIONS 

We made computer simulations of the model. The results accorded with the theoretical 
predictions quite well at most parameters ranges. Figure ^a), (b) and (c) give simulations 
of networks with size N = 1000000, m = 5 and / = —2, 2, 5 respectively. With log-log 
scale, the main part of the distributions give slope lines in the figure. The slopes of linear fit 
approximately equal to the theoretical results (7 = 2.6, 3.4, 4 respectively). The prediction of 
the ratios of x degree to y degree are also confirmed by numerical simulations (see figure |2J) . 
One should notice that in the figures, we had made shifts for x degree distributions and total 
degree distributions from P x (k) and P x+y (k) to P x (k — f) and P x+y (k — f) respectively. 
Therefore, their log-log plots behave as straight lines parallel to those of y degree. 

However, we found that when we take small g, the simulations are always different from 
theoretical predictions. With small g, when / is positive, the component of x degree is 
larger than theoretical value while the component y degree accords with theoretical value 
(see figureEfa)); and when / is negative, the component oiy degree is larger than theoretical 
value while the component of x degree accords with theoretical value(see figure Efb)). In 
theoretical calculus, at infinite t limitation, we neglect the last term of equation JSJ). When 
g is small, decay of this term is slow, and actually in numerical simulations, t is not infinite. 
Therefore, the neglect of this term make the bias of simulation results from predictions at 
small g. 

We also looked at the fluctuation of x degree Xk to theoretical value as a function of total 
degree k. It illustrates how the competition generates heterogeneity in edge composition of 
each vertex. We calculated the relative standard deviation respectively to theoretical value. 
The fluctuation decreases quickly with the increase of total degree nearly in a power law 
way (see figure EJ). We did not calculate the standard deviation of Xk/yu because y degree 
may be zero for some vertexes. 

IV. DISCUSSIONS 

The mathematical expressions of degree distribution of both Y sub-network and X sub- 
network provide deep insight into the dynamics of evolving systems. We build a competitive 
environment where not only nodes but also different types of links compete. This model 
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can reflect many properties of social network. A newly added node has a fixed number 
of m adding edges. However, other nodes decided how many X edges and Y edges there 
would be. Obviously, other nodes compete for these m edges by what they already have. 
For instance, let X-edge denotes financial relationship between individual persons, and Y 
represents other connection. Rich people tend to have more financial relationship with other 
people, while interestingly, more financial links signifies his richness. But when we focus on 
other links between people, rich people are not necessarily to be so lucky. Due to the 
limitation of personal capability, time, devotion, one can not have infinite connections with 
others. Therefore, it is at the expense of less Y edges to get more X edges, and vice versa. 
The above study shows the relationship between X and Y edges. 

An important character of this model is that the degree distribution of X sub-network 
shows a linear shift while Y sub- network does not. Simulations confirm this character, 
though the exact value of the shift may has some errors, due to the application of mean-field 
approximation. The difference of linear shift between X and Y sub-network comes from 
different preferential attachment. 

In the following, we will discuss how the variables g, f and m work on in our theory. 

m is always thought as a fixed number in a certain network, for we can identify at least 
vaguely how many links are added in each time interval. The value of / varies from — m to 
oo, but if / is too large both of total degree distribution and x degree distribution approach 
exponential increment, 

(jfc + f)'^+i) = e -(3+£)ln(fe+/) = e -(3+£)pn/+ln(l+$)] 

w f-( 3 +^e^ 3+ ^7 = /-( 3 +£) e -^ + 7) (17) 

And if / >> m, the characteristic degree of total network and X sub- network is m. On the 
other hand, when / is taken a bit large, the actual probability to connect Y links will be so 
small that in the end most of the nodes have small y degree, while the rest have very large 
y degree. 

/ rng \ 

Note that in the several derivations above, we neglect the term Const -t l2 ( 2m +fl' for time 
t is considered approaching to infinite. However, this ideal condition can not be reached to 
apply the model, i.e. t will always be finite actually, even considerable large. So problem 
rises that if g, f and m are chosen or set so that the component 2 (2m+f) ' ls sman enough and 
the whole term becomes an unneglectable factor to the degree distribution. The derivation 
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above may not hold true. But if we take this term into account, we will find ourself awkward 
in searching a precise solution for the function of degree distribution. 

The most intriguing parameter in the parameters is g. Apparently, from above discussion, 
g should not be chosen too small to make the term t v W^+7)> unneglectable. The simulations 
showed that small g make the degree distributions depart from predictions, g should not be 
too large, as well. A very large g in finite time lead the network to have few Y links also, as 
we have found by computer simulation. Based on above reasons, g is preferred to be such a 
number that is neither too large, nor too small. Here we only use vague words "large" and 
"small", for these limitations are due to the finity of time a network evolves and the finity 
of nodes and links it contains. To say a bit more accurate, we need g to be big enough, with 

a view to t, to make the term t y W^+n> small enough to be neglected; but g should not be 

, , , , dxi Xi + yi + f + gy { dy { x { + y { + f - gy { 

too larger to make — — = m —r too larger than — — = m 



dt 2Q> i + %• + /) ° dt 2 (J2x j + y j + fy 

which leads the network having few Y edges. 

For a given node with degree k at given time, the probability p x (k) that how many X 
edges it has follows equation (JH). The figure |U shows the fluctuation decreases quickly with 
k. It seems that the distribution of fluctuation should follows the binomial distribution. 
However, here the degree of a node is varying with time. It increased from m to k and p x 
varied synchronously. It results in the fact that the fluctuation of the number of X edges is 
much smaller than that of the binomial distribution. 

Finally we discuss directed network and find some difference between this kind of network 
and what we have studied. We give general results below. Every time step, we introduce m 
directed edges. The whole degree has a m increment (see j^)- 

J2(h + f) = (m + f)t 

using the same approaches, we get 

P x+y (k) = (1 + {k + (lg) 

And the degree distributions of X subnetwork and Y subnetwork are 

PM = (i + l ) [ m tl ] ~ l ^ +i) ^ + fr {2+ - ] ( 19 ) 

V mj 2+g 
P y (k) = (l + £) (-L)^)*-^) (20) 
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The ratios of total degree to X degree and Y degree are 



P x+y (k-f):P x (k-f):P y (k) 
/ /(2 + g)-l 




f(2 + g) 



1 




Compared with the results of undirected network (JTfjJ), one can find that proportion of Y 
edges of directed network is larger than that of undirected network. 
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FIG. 1: The distribution of total degree, X degree and Y degree. The size of the network N = 
1000000 and m = 5. The lines are linear fits of the main part of the data, (a) / = —2, g = 3. 
The slopes of the lines are j k = —2.5, 7 X = —2.5 and "y y = —2.6. The prediction is 7 = —2.6. (b) 
/ = 2, g = 1. The slopes of the lines are 7 fc = —3.3, 7 X = —3.4 and 7 y = —3.4. The prediction is 
7 = —3.4. (a) / = 5, g = 2. The slopes of the lines are 7 fc = —3.8, r y x = —4.0 and 7^ = —3.9. The 
prediction is 7 = —4. 
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FIG. 2: The ratios of the distribution X degree to the distribution of Y degree with respective 
to different g. The curves are theoretical predictions. The size of the network N = 1000000 and 
m = 5. (a) / = -2; (b) / = 2; (c) / = 5. 
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FIG. 3: The degree distribution of total degree, x degree and y degree with small g = 0.1. The 
size of the network N = 1000000 and m = 5. (a) / = -2. (b) / = 2. 
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FIG. 4: The relative standard deviation of X degree respective to total degree. The size of the 
network N = 1000000. 
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